a=[2,-3,1;8,3,2;45,1,-9]; 
b=[4;2;17]; 
x=inv(a)*b
